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The new EGRET observations of PSR B1706-44 have not yet been successfully 
exploited to discriminate among the various proposed models which attempt to 
describe the location and mechanisms for 7 -ray emission from pulsars. They have 
stimulated (still unsuccessful) attempts to fit these data, together with those from 
other 7 -ray pulsars, (cf., §5.5 of D.J. Thompson et al. in the enclosed reprint - the 
collaboration’s publication), and to understand the absence of strong pulsed 7 -ray 
sources among the millisecond pulsar candidates. 

A particularly puzzling feature is the success of a single power law extrapola- 
tion (dN/du oc u~ 1 - 7 ) photon intensity from GeV 7 -rays through many orders of 
magnitude, in some cases all the way down to the 10 2 MHz - 1 GHz microwave 
regime (except, of course, for Geminga where no radioemission has been detected). 
PSR 1706 is remarkable in that this power law fit extends to above 10 GeV. It is 
suggestive of 7 -ray production mechanisms quite different from those which have 
been the basis of most 7 -ray pulsar models. If so, the present data from PSR B1706 
and its siblings may at last reveal how and why pulsar accelertors are formed. Un- 
fortunately, the answer is not so obvious that this is now the case. A reprint of 
the collaboration’s published work is here attached (D.J. Thompson et al. 1996). 
Clearly the observations are far more complete and robust than are the interpreta- 


tions of them. 
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ABSTRACT 

The Energetic Gamma-Ray Experiment Telescope (EGRET) on the Compton Gamma Ray Observatory 
has observed PSR B1706 — 44 a number of times between 1991 and 1995. From these data, a more 
detailed picture of the gamma radiation from this source has been developed, showing several character- 
istics that distinguish this pulsar from others: the light curve is complex, with evidence of three pulses; 
there is no detectable unpulsed emission; and the energy spectrum changes slope at 1 GeV. No two of 
the known gamma-ray pulsars have quite the same characteristics; this diversity makes interpretation in 
terms of theoretical models difficult. 

Subject headings : gamma rays: observations — pulsars: individual (PSR B 1706 — 44) 


1. INTRODUCTION 

In 1981, the second COS B catalog of high-energy 
gamma-ray sources included an entry for 2CG 342 — 02, an 
unidentified E > 100 MeV source near the Galactic coordi- 
nates indicated by the source name (Swanenburg et al. 
1981). Nearly 10 years later, radio pulsar PSR B1706— 44 
was discovered during a high-frequency search of the 
southern sky (Johnston et al. 1992). These lines of research 
from opposite ends of the electromagnetic spectrum con- 
verged shortly thereafter with the identification, using data 
from the Energetic Gamma-Ray Experiment Telescope 
(EGRET) on the Compton Gamma Ray Observatory 
(CGRO), that the high-energy gamma radiation was pulsed 
at the radio period (Thompson et al. 1992). 

EGRET has now detected pulsed high-energy gamma 
radiation from six pulsars (PSR B 1706 — 44 and Crab, 
Nolan et al. 1993; Vela, Kanbach et al. 1994; PSR 
B1951 + 32, Ramanamurthy et al. 1995; Geminga, Bertsch 
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et al. 1992, Mayer-Hasselwander et al. 1994; and PSR 
B1055 — 52, Fierro et al. 1993), with some evidence for a 
seventh, PSR B0656 + 14 (Ramanamurthy et al. 1996). 
Upper limits have been calculated for selected samples of 
radio pulsars (Thompson et al. 1994; Fierro et al. 1995) and 
for all cataloged pulsars (Nel et al. 1996). In particular, PSR 
B 1509 — 58, which is seen at energies up to 1 MeV (Wilson 
et al. 1993 ; Matz et al. 1994; Carraminana et al. 1995), is not 
detected by EGRET (Brazier et al. 1994). The present work 
is a detailed analysis of the EGRET observations of PSR 
B1706 — 44, based on repeated observations during 1991— 
1995 that have tripled the source exposure time compared 
with the discovery data. 

The basic pulsar parameters (Taylor, Manchester, & 
Lyne 1993) are shown in Table 1. McAdam, Osborne, & 
Parkinson (1993) have suggested an association of PSR 
B 1706 — 44 with a shell-type supernova remnant (SNR), 
G343.1 — 2.3. Based on distance estimates that may conflict, 
and on the lack of interaction of the pulsar with the 
remnant, Frail, Goss, & Whiteoak (1996) consider the 
association unlikely, though possible. They suggest instead 
that the pulsar lies within a small plerion seen in their 20 cm 
radio image. Unpulsed X-ray (Becker, Brazier, & Triimper 
1995) and TeV (Kifune et al. 1995) emission have been 
detected from this source, although pulsations have not 
been reported at either of these wavelengths. As noted by 
Becker et al. (1995), the pulsed X-ray upper limit of 18% 
from ROSAT does not rule out a pulsed fraction similar to 
other pulsars observed by ROSAT (e.g., Vela is 11% pulsed; 
Ogelman, Finley, & Zimmerman 1993). 

2. RADIO TIMING OBSERVATIONS 

Shortly after its discovery, PSR B 1706 — 44 was added to 
the list of nearly 300 pulsars that were monitored regularly 
by radio astronomers to assist gamma-ray telescopes on the 
CGRO (Johnston et al. 1995). High-energy gamma-ray data 
are sparse; weak, short-period gamma-ray pulsars are only 
detectable if the timing parameters are determined indepen- 
dently of the gamma-ray data. In the case of PSR 
B 1706 -44, this monitoring, carried out at Parkes Obser- 
vatory in New South Wales, has continued. A glitch with 
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TABLE 1 
Basic Data 


Parameter Value 


Names PSR B1706 - 44, PSR J1709 - 44, PSR J1709-4428 

Period P 0.1024 s 

Period derivative P 93.0 x 10“ 15 s s -i 

Timing age t 17,000 yr 

Spin-down luminosity is 3.4 x 10 36 ergs s -1 

Magnetic field B 3.1 x 10 12 G 


<5 P/P ~ 2 x 10“ 6 occurred at MJD = 48775 + 15 
(Johnston et al. 1995). In addition to this glitch, the pulsar 
exhibits considerable timing noise. For this reason, the 
timing solutions used for the EGRET analysis were devel- 
oped piecewise over time intervals for which the pulse phase 
could be adequately modeled using only a simple spin- 
down law. Table 2 lists the solutions relevant to the CGRO 
viewings, given in terms of frequency v and its derivatives 
instead of period, and valid at time T 0 . These timing solu- 
tions are from the database maintained at Princeton Uni- 
versity (anonymous FTP: ftp://puppsr.princeton.edu/gro). 

Figure 1 shows the difference between the fit and the 
measured radio arrival time (timing residuals) for one of 
these solutions. The fact that the observed radio pulse never 
varies from the solution by even as much as 0.2 ms is a 
strong indication that the gamma-ray arrival times can be 
converted to pulsar phase with comparable accuracy. 
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Fig. 1. — Timing residuals for PSR B 1706 — 44. The residual is the dif- 
ference between the modeled and observed radio pulse arrival time. 


3. GAMMA-RAY OBSERVATIONS 

EGRET, the high-energy gamma-ray telescope on the 
CGRO , is described by Thompson et al. (1993, and refer- 
ences therein). Operating from about 30 MeV to over 20 
GeV, EGRET records individual gamma-ray photons as 
electron-positron pair production events. Each photon is 
analyzed to provide its arrival direction, its energy, and its 
arrival time (in Universal Coordinated Time [UTC]) to an 
accuracy of better than 100 ps. The field of view mapped by 
EGRET extends to more than 30° from the instrument axis, 
although the sensitivity at angles beyond 30° is less than 
15% of the on-axis sensitivity. Because of the low flux level 
of the high-energy gamma rays, viewing intervals are typi- 
cally 2-3 weeks. 

PSR B 1706— 44 was within 30° of the telescope axis 
during 23 of the CGRO viewing intervals to date. No addi- 
tional EGRET observations of this pulsar are scheduled. 
Table 3 is a list of the EGRET observations, showing the 
dates, the angle between the pulsar and the EGRET axis, 
and the EGRET exposure for the energy range greater than 
100 MeV. 

Although PSR B1706 — 44 is a relatively bright source 
compared with many others seen by EGRET (cf. the second 
EGRET catalog; Thompson et al. 1995), its gamma-ray 
count rate is nevertheless low, about two photons (E > 100 
MeV) per hour of exposure when the source is within 10° of 
the EGRET axis. Data processing for EGRET relies on two 
principal methods: timing analysis and spatial analysis. 

For the timing analysis, a set of detected photons is selec- 
ted within an energy-dependent cone of radius 

6 < 5?85(£ y /100 MeV) -0 534 , (1) 

with respect to the pulsar position (E y in MeV). This choice 
represents the energy-dependent angular resolution of the 
EGRET instrument (Thompson et al. 1993). The conversion 
of gamma-ray arrival time at the location of the CGRO to 
pulsar phase is carried out using a modification of the 
TEMPO timing program (Taylor & Weisberg 1989) and the 
JPL DE200 ephemeris. 


TABLE 2 

Radio Timing Parameters for PSR B1706-44 


T 0 v v v 

Valid Dates (MJD) (s* 1 ) (1(T 12 s“ 2 ) (10 22 s -3 ) 


1990 Jan 18-1992 May 3 48327 9.7612721442486 -8.86376 1.21 

1992 Jul 9-Dec 3 48885 9.7608645998410 -8.89014 16.6 

1993 Jan 13-Jul 8 49088 9.7607088316970 -8.87504 6.08 

1993 Jul 19-1994 Dec 22 49447 9.7604336985587 -8.86669 2.19 

1994 Aug 25-1995 Jul 22 49754 9.7601985856255 -8.86125 1.97 
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TABLE 3 


EGRET Observations of PSR B 1706 -44 


Viewing Number 

Observation Dates 

Aspect 

Angle 

Exposure 
(>100 MeV) 
(10 7 cm 2 s) 

0050 

1991 Jul 12-26 

16?9 

25.9 

0160 

1991 Dec 12-27 

28.3 

10.8 

0230 

1992 Mar 19-Apr 2 

21.7 

5.7 

0270 

1992 Apr 28-May 7 

12.0 

9.7 

0350 

1992 Aug 6-11 

24.1 

1.9 

0380 

1992 Aug 27-Sep 1 

24.1 

3.0 

2100 

1993 Feb 22-25 

15.4 

2.3 

2140 

1993 Mar 29-Apr 1 

15.4 

3.1 

2190 

1993 May 5-8 

19.8 

0.9 

2230 

1993 May 31-Jun 3 

16.2 

3.0 

2260 

1993 Jun 19-29 

14.2 

10.5 

2290 

1993 Aug 10-11 

23.2 

0.5 

2295 

1993 Aug 12-17 

23.2 

2.6 

2320 

1993 Aug 24-Sep 7 

5.2 

19.0 

3023 

1993 Sep 9-21 

21.8 

5.8 

3230 

1994 Mar 22-Apr 5 

16.1 

16.2 

3340 

1994 Jul 18-25 

26.4 

3.3 

3365 

1994 Aug 4-9 

6.2 

7.7 

4210 

1995 Jun 6-13 

12.6 

8.8 

4220 

1995 Jun 13-20 

12.5 

11.2 

4230 

1995 Jun 20-30 

19.6 

10.6 

4235 

1995 Jun 30- Jul 10 

16.4 

11.3 


The spatial analysis compares the observed gamma-ray 
map with that expected from a model of the diffuse Galactic 
and extragalactic radiation (Hunter et al. 1996). Source 
location and flux as a function of energy are determined 
using a maximum likelihood method (Mattox et al. 1996). 
The timing and spatial approaches can also be combined to 
produce phase-resolved maps and energy spectra. 

4. RESULTS 
4.1. Light Curve 

Compared with the very bright Crab, Vela, and Geminga 
gamma-ray pulsars (Ramanamurthy et al. 1995), PSR 
B1706 — 44 is considerably weaker and lies in the Galactic 
center region where the diffuse Galactic gamma radiation is 
much more intense. For these reasons, the results are char- 
acterized by somewhat larger uncertainties than seen for 
these other pulsars. 

Figure 2 shows the light curve for PSR B 1706 — 44 for 
energies above 400 MeV, including all viewing periods 
listed in Table 3. The peak of the single radio pulse 



Pulsar Phase 


Fig. 2. — High-energy gamma-ray light curve for PSR B1706 — 44. The 
102 ms period is divided into 50 phase bins. The single radio pulse defines 
the zero phase. 


(Johnston et al. 1992) defines phase 0 on this 50 bin plot. 
The gamma-ray energy selection is based on maximizing 
the significance of the pulsed signal, as characterized by the 
(binned) x 2 or (unbinned) H-test value (de Jager, Swanepoel, 
& Raubenheimer 1978). For lower energies, the larger 
acceptance cone introduces more Galactic background, 
while the statistics are reduced for higher energies. 

The gamma-ray light curve, which for this selection con- 
tains 1325 gamma rays, is more complex than those of the 
Crab, Vela, Geminga, and PSR B 195 1 H- 32, all of which 
show two narrow peaks separated by 0.4-0. 5 in phase. In 
the case of Vela, for example, the pulse widths (FWHM) are 
a few milliseconds (Kanbach et al. 1994). PSR B1055 — 52 
(Fierro et al. 1993) has what appears to be a single broad 
pulse in gamma rays, but the statistics for this pulsar are far 
smaller than those for PSR B 1706 — 44. Between phases 0.2 
and 0.6, PSR B 1706 — 44 shows at least two, and possibly 
three, peaks. Figure 3 shows the same data plotted in 100 
bins (and explicitly displaying the uncertainties due to 
counting statistics) with three possible fits to the pulse 
shape: (1) a fit consisting of a single Gaussian peak plus a 
constant background (x 2 = 139 for 96 degrees of freedom 
[dof] [100 bins minus 4 degrees of constraint for the fit], 
with probability = 0.007 for this x 2 ); (2) a fit using two 
Gaussians plus a constant (x 2 = 99 for 93 dof, prob- 
ability = 0.32); and (3) a fit with three Gaussians plus a 
constant (x 2 = 90 for 90 dof, probability = 0.48). Although 
the fit using a single peak can be rejected with high con- 
fidence, the difference between the two Gaussian and three 
Gaussian fit is not compelling. An F-test indicates that the 
two Gaussian plus constant fit has a 60% probability of 
being acceptable. Nevertheless, the third (middle) peak 
appears to be a persistent feature, appearing consistently 
when the data set is divided in half by time (Fierro 1995). 
Based on these fits, the two principal peaks in the light 
curve fall at phases 0.27 + 0.02 and 0.45 ± 0.02 and have 
FWHM of approximately 12 ms each. The third peak, if 
included, falls at phase 0.36 + 0.02 and has FWHM = 3 ms. 

Figure 4 shows the light curve in four energy bands. The 
decrease in background and the decrease in statistics with 
increasing energy are both clear from this figure. Above 5 
GeV, all 12 photons fall into the pulse region. The basic 



Fig. 3. — High-energy gamma-ray light curve for PSR B1706-44, with 
100 phase bins. Uncertainties due to counting statistics are shown. The 
fitted curves, as described in the text, are two Gaussian peaks plus a 
constant background ( heavy solid line), three Gaussian peaks plus a con- 
stant background {light solid line), and single Gaussian peak plus a con- 
stant background ( heavy dotted line). 
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Pulsar Phase 


Fig. 4. — Gamma-ray light curve for PSR B 1706 — 44 in four energy 
ranges, with 50 bin resolution. 


structure of the light curve is similar at all energies, with 
strong indications of two peaks and some evidence of three 
peaks seen at all energies. The evidence of pulsed emission 
extends well above 5 GeV, in contrast to Geminga and 
Vela, which show strong cutoffs above a few GeV. Of the 12 
highest energy photons, five have measured energies above 
10 GeV, and one has a measured energy of 64 ± 34 GeV. 

From 30 to 100 MeV, the pulsed signal is harder to 
separate from the surrounding Galactic background. The 
signal is still visible, however, if a narrower cone is taken 
rather than the one described above. The EGRET point- 
spread function itself has a narrow and a broad component 
(Thompson et al. 1993), and by selecting a narrower range, 
specifically a fixed cone of half-angle 2? 5, it is possible to 
obtain a light curve whose probability of occurring by 
chance is less than 10 ~ 5 . The pulse shape is consistent with 
that seen at the higher energies. 


4.2. Search for Unpulsed Emission and Time History of 
T otal Emission 

Although Figure 4 makes it clear that the gamma radi- 
ation from PSR B 1706 — 44 is strongly dominated by the 
pulsed emission at high energies, a search for unpulsed 
emission can be undertaken over broader energy ranges. In 
the region of sky mapped by EGRET around the pulsar, 
each photon's arrival time is converted to pulsar phase, 
whether or not this photon is likely to have come from the 
pulsar itself. Phase-resolved maps of the sky are then con- 
structed, and the spatial analysis using maximum likelihood 
(Mattox et al. 1996) is used to assess the statistical signifi- 
cance of a source at the pulsar location. The likelihood ratio 
test is used to determine the significance of point sources. 
The likelihood ratio test statistic is TS = 2(lnLi — lnL 0 ), 
where In L v is the log of the likelihood of the data if a point 
source is included in the model, and In L 0 is the log of the 
likelihood of the data without a point source. For positive 
values of TS, the analysis gives the most likely gamma-ray 
flux of a source at the pulsar location. The result in each 
phase bin for E > 100 MeV is shown in Figure 5. The pulsar 
is detected with high significance only between phases 0.2 
and 0.55, exactly the range expected from the light curves. 
There is weaker evidence for a leading wing and a trailing 
wing, but no evidence for a source in the rest of the phases. 
Table 4 shows the results of the phase-resolved analysis for 
two broad energy ranges. In each case, the flux is time 
averaged over the entire phase. In the absence of any 
unpulsed emission, the pulsed flux is equivalent to the total 
flux. Above 100 MeV, the unpulsed component is less than 
20% of the pulsed (total) emission. 

Figure 6 shows the E > 100 MeV EGRET observations 
of PSR B1706— 44 flux as a function of time, from 1991 to 



Fig. 5 . — E > 100 MeV gamma-ray flux at the location of PSR 
B 1706 —44 as a function of pulsar phase. 


TABLE 4 

Pulsed and Unpulsed Results 


Unpulsed Limit 

Energy Range Pulsed Flux (95% confidence) 


30 < £ < 100 MeV 1.05 ± 0.24 x 10“ 6 <1.01 xl0“ 6 

£>100 MeV 1.28 ± 0.08 x 10“ 6 <0.24 x 10“ s 
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Fig. 6 . — E > 100 MeV gamma-ray flux at the location of PSR 
B1706— 44 as a function of time. Each data point represents one viewing 
period. The uncertainties shown are statistical only. 



0 4 0.6 

Pulsar Phase (Radio pulse = 0.0) 


Fig. 8 —Hardness ratio (flux for E > 300 MeV/flux for 100-300 MeV) 
of PSR B 1706 — 44 as a function of pulsar phase. 


PSR B 1706 — 44 appears to be steady. The x 2 is 13 for 13 
dof. 


1994. The 1995 observations were affected by known drifts 
in the EGRET performance due to the deterioration of the 
spark chamber gas and have therefore been excluded. The 
dashed line is the flux calculated from the summed exposure 
to the pulsar. As for the Crab, Geminga, and Vela pulsars 
(Ramanamurthy et al. 1995), the gamma radiation from 



Fig. 7. — Phase-averaged energy spectrum of PSR B 1706 — 44. The fit to 
a broken power law is described in the text. The uncertainties shown are 
statistical only. 


4.3. Energy Spectrum 

Because there is no evidence for unpulsed emission, the 
energy spectrum of the pulsar can be derived by defining the 
background level in any energy band as the flux seen in the 
off-pulse part of the light curve. The excesses in each band 
are then compared with model spectra forward folded 
through the EGRET energy response function, as described 
by Nolan et al. (1993). Pulsed emission is detected from 50 
MeV to more than 10 GeV. The spectrum, shown in Figure 
7 as a phase-averaged photon number spectrum, is well 
described by a broken power law, 

dN ( E V 

dE - 9 \imzvj 

photons cm -2 s" 1 MeV' 1 , (2) 

where 

f- 1.27 ± 0.09E < 1000 MeV , 

* ~ ( — 2.25 + 0.13E ^ 1000 MeV . (3) 

The Crab, Vela, and Geminga pulsars all show spectral 
variations as a function of pulsar phase in the EGRET 
energy range. Although PSR B 1706 -44 is weaker than 
these three, the possibility of spectral variation can be 
examined. Figure 8 shows the hardness ratio, the flux above 
300 MeV divided by the flux from 100 to 300 MeV, for PSR 
B 1706 — 44 as a function of phase, using the same 20 phase 
bins as in Figure 4. With the possible exception of the 
leading wing of the pulse (phase 0.15-0.20), the hardness 
ratio is consistent for all the pulsed emission. Energy 
spectra constructed for portions of the phase corresponding 
to the peaks are also consistent with each other (Fierro 
1995). 

5. DISCUSSION 
5.1. Light Curves 

The overall gamma-ray light curve for PSR B1706 — 44 
differs markedly from those of most of the other gamma-ray 
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pulsars. The Crab, Vela, Geminga, and PSR B 195 1 + 32 
light curves are all chracterized by two narrow pulses 
separated by 0.4-0.5 in phase. PSR B 1509 -58 (seen only 
< 1 MeV) has a well-defined single pulse, and PSR 
B 1055 — 52 has a broad single pulse (although not well 
defined due to limited statistics). PSR B 1706 — 44 shows two 
broader pulses with a phase separation of 0.2. This distinc- 
tion points strongly to a different beaming geometry for this 
pulsar than for the others. Radio polarization studies, for 
example, indicate that PSR B 1706— 44 has its rotation and 
magnetic axes more aligned than Vela (Guojun et al. 1995). 
The possible third peak in the light curve is unique to PSR 
B1706 — 44 among gamma-ray pulsars. In this respect, it 
resembles radio pulsar light curves that have both core and 
conal emissions (e.g., Rankin 1991). 

5.2. Distance 

The distance determined from the dispersion measure 
(DM) 75.69 cm' 3 pc (Taylor & Cordes 1993) is 1.8 ( — 0.5, 
+ 0.7) kpc. Kinematic distance limits derived from H i 
absorption (Koribalski et al. 1995) are D Lowcr = 2.4 + 0.6 
kpc and Z> Upper = 3.2 + 0.4 kpc. The possible association of 
PSR B 1706 -44 with the SNR G343.1 -2.3 (McAdam et al. 
1993) would imply a distance of 3 kpc based on the X-D 
relation. For purposes of calculation, we adopt a distance of 
2.4 kpc, consistent with both the DM and H i measure- 
ments. It is important to recognize that all such distance 
estimates have uncertainties of 25% or more. 

5.3. Luminosity and Beaming 

In terms of the observed energy flux F £ , the gamma-ray 
luminosity of a pulsar is 

L y = 4nfF E D\ (4) 

where / is the fraction of the sky into which the pulsar 
radiates and D is the distance to the pulsar. This beaming 
fraction / is highly uncertain. Based on the pulse shape, 
Thompson et al. (1992) used a value of 0.27. In a nearly 
aligned rotator model, Dermer & Sturner (1994) find a 
beaming fraction of less than 0.1, while an outer-gap model 
(Yadigaroglu & Romani 1995) suggests a value of 0.55. In 
comparing the EGRET-detected pulsars, Thompson et al. 
(1994) adopted a value of 1/4 tl 

The observed energy flux obtained by integrating equa- 
tion (2) in the range 100 MeV- 10 GeV is (9.4 + 1.2) 
x 10” 10 ergs cm -2 s. For a distance of 2.4 kpc, the gamma- 
ray luminosity of PSR B 1706 -44 is then (5 + 3) 
x 10 34 x 4nf ergs s -1 . Unless the beaming fraction is ex- 
tremely small, the gamma radiation represents on the order 
of 1% of the spin-down luminosity, E = 3.4 x 10 36 ergss _A . 

5.4. Broadband Energy Spectrum 

The multiwavelength spectrum of PSR B 1706— 44 is 
shown in Figure 9, using two different formats. At the top, 
the photon number spectrum is shown. The photon spec- 
trum falls steadily across the entire electromagnetic spec- 
trum. At the bottom, the frequency versus vF v display 
shows the observed power per log (frequency). Pulsed emis- 
sion is detected from PSR B 1706 -44 only in the radio and 
high-energy gamma-ray bands. Unpulsed emission, which is 
associated with the pulsar on the basis of positional infor- 


Energy (keV) 



10 9 10' 2 10 15 10 ,B 10 21 10 24 10 27 
Frequency (Hz) 

Fig. 9. — Multi wavelength spectrum of PSR B1706 - 44. Top , photon 
number spectrum; bottom, the frequency v vs. vF v . References are radio 
(Taylor et al. 1993; Johnston et al. 1992); X-ray ( ROSAT ), the 2 a pulsed 
upper limit is 18% of the unpulsed emission (Becker et al. 1995); gamma 
ray, OSSE (Schroeder et al. 1995) and COMPTEL (Carraminana et al. 
1995); EGRET (present results); TeV detection (unpulsed) (Kifune et al. 
1995); and TeV pulsed upper limit (Nel et al. 1994). 


mation, is seen in the ROSAT X-ray range (Becker et al. 
1995) and in the very high energy gamma-ray regime using 
ground-based detection of the Cerenkov radiation from 
gamma-ray-induced particles (Kifune et al. 1995). Because 
of the very flat energy spectrum below 1 GeV, the OSSE 
(Schroeder et al. 1995) and COMPTEL (Carraminana et al. 
1995) observations yield only upper limits. EGRET obtains 
positive detections from 50 MeV to about 20 GeV. The 
maximum vF v occurs near 1 GeV, the break in the observed 
EGRET spectrum. The observable power from this pulsar is 
strongly dominated by the gamma-ray emission. 

5.5. Gamma-Ray Pulsar Models 

Two general classes of models have been proposed for 
gamma-ray pulsars. In polar cap models (recent examples 
are Daugherty & Harding 1994, 1996 and Sturner & 
Dermer 1995), the particle acceleration and gamma-ray 
production takes place in the open field line region above 
the magnetic pole of the neutron star. In outer gap models 
(recent examples are Ho 1993 and Romani & Yadigaroglu 
1995), the interaction region lies in the outer magnetosphere 
in vacuum gaps associated with the last open field line. In 
either of these models, PSR B 1 706 — 44 could be expected to 
be a detectable gamma-ray pulsar, because its physical 
parameters are quite similar to those of Vela. The details of 
the gamma-ray observations, however, are not all readily 
explained by either model. 

Because either type of model can be viewed as having a 
hollow cone geometry, a double pulse has a straightforward 
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explanation. The observer’s line of sight cuts across the edge 
of the cone at two places. Although the specific details 
depend on the size of the beam, and its relationship to the 
rotation axis and the observer’s line of sight, in the case of 
PSR B 1706 — 44, one possibility is that the line of sight is 
closer to the edge of the cone than for the pulsars with two 
widely spaced light-curve peaks. The fact that the peaks are 
broader for PSR B 1706 — 44 is also consistent with this geo- 
metric picture, because the line of sight crosses the cone at a 
shallower angle. 

Following from the geometric picture, the existence of a 
third pulse is not expected in the simplest version of either 
model. An earlier version of the outer gap model (Cheng, 
Ho, & Ruderman 1986) predicted that gamma-ray pulsars 
could have up to four peaks, with the second pair of peaks 
arising from particle beams moving inward as well as 
outward along the outer gap. More recent outer gap models 
expect much less radiation from the inward-moving beams 
(Chiang & Romani 1994). A unified polar cap/outer gap 
model (Kamae & Sekimoto 1995) can explain three peaks in 
the light curve, but this model predicts that the middle peak 
should be broad and have a soft spectrum, in contrast to the 
present observations. In the absence of any evidence for a 
third pulse in the other gamma-ray pulsars, PSR B 1706 — 44 
appears to be distinctive in this respect. 

The energy spectrum presents similar challenges. In the 
polar cap models, a sharp turnover is expected in the few to 
10 GeV energy range due to attenuation of the gamma-ray 
flux in the magnetic field (Daugherty & Harding 1994, 
1996). The spectral change in the 1 GeV range does not 
appear to fit this model, although there could be a sharp 
break in the energy range above 20 GeV, as suggested by 
the absence of pulsed TeV emission. On the other hand, the 
lack of variation of the spectrum with pulsar phase for PSR 
B1706 — 44 is expected in the Daugherty & Harding (1996) 
polar cap model. The outer gap model does predict a flat- 
tening of pulsar spectra at lower energies (Ho 1993), 
although, in the case of Vela, this flattening occurs well 
below 100 MeV, not near 1 GeV, and the flattening is to an 
E ~ 0 61 spectrum, not E~ 1 ’ 27 as observed. The dependence 
of the spectrum on the geometry has not been extensively 
studied in the outer gap model. It is fair to say that the 
spectral shape of PSR B 1706 — 44 was not predicted by any 
model. 

A reasonable hypothesis is that the difference in the 
beaming/observer geometry (compared with the brighter 
gamma-ray pulsars), which produces the unusual pulse 
shape for PSR B 1706 — 44, is also responsible for the differ- 
ent energy spectrum. How these observed features relate to 
the gamma-ray production process is unclear. PSR 
B1055 — 52 is the pulsar closest to PSR B1706 — 44 in light- 
curve shape; therefore, its energy spectrum might show a 
similar feature. Unfortunately, the statistics for both the 
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light curve and the energy spectrum of PSR B1055 - 52 are 
much more limited (Fierro et al. 1993) and cannot at present 
address this question. 

6. SUMMARY 

PSR B1706 — 44 is one of at least six spin-powered pul- 
sars seen at gamma-ray energies above 100 MeV. Obser- 
vations from EGRET on the CGRO between 1991 and 
1995 have provided new details of the gamma radiation : 

1. The light curve is complex, with at least two peaks 
separated by about 0.2 in phase, with some evidence of a 
third peak. None of the other gamma-ray pulsars show a 
similar light curve. 

2. There is no detectable unpulsed gamma radiation 
from the pulsar. 

3. There is no evidence that the gamma-ray flux from the 
pulsar varies on long timescales. 

4. The gamma-ray energy spectrum can be represented 
by a broken power law, with a break near 1 GeV. The 
spectral index changes by 1.0, in contrast to the Vela and 
Geminga spectra, which have sharp spectral breaks (index 
change greater than 1) in the few GeV energy range. 

5. The maximum observable power from the pulsar is in 
the gamma-ray energy range. The observed X-ray emission 
is unpulsed, with a pulsed fraction less than 18%. The TeV 
emission is also not pulsed. The Crab is the only other 
pulsar seen at TeV energies, and its TeV emission is not 
pulsed. 

6. The gamma radiation represents about 1% of the spin- 
down luminosity of the pulsar, although the unknown 
beaming geometry makes this estimate uncertain. 

7. Present gamma-ray pulsar models offer no straightfor- 
ward explanation of either the light-curve shape or the 
energy spectrum. 
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